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Observation of near-critical reflection of internal waves 
in a stably stratified fluid 
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An experimental study is reported of the near-critical reflection of internal gravity waves over 
sloping topography in a stratified fluid. An overturning instability close to the slope and triggering 
the boundary-mixing process is observed and characterized. These observations are found in good 
agreement with a recent nonlinear theory. 

PACS numbers: 47.55.Hd,47.35.+i,47.20.-k 



I. INTRODUCTION 

The reflection of near-critical internal waves over slop- 
ing topography plays a crucial role in determining ex- 
changes between the coastal ocean and the adjacent deep 
waters. Direct measurements of mixing in the ocean, us- 
ing tracers Q, have vindicated decades of phenomeno- 
logical and theoretical inferences. In particular, these 
measurements have shown that most of the vertical mix- 
ing is not taking place inside the ocean, but close to the 
boundaries and topographic features Q. These results 
have directed attention to the possible role of internal 
wave reflection in the boundary-mixing process. 

Internal waves have different properties of reflection 
from a rigid boundary than do sound or light waves Q. 
Instead of following the familiar Snell's law, internal 
waves reflect off a boundary such that the angle with 
respect to gravity direction is preserved upon reflection 
(Fig. GJ. This peculiar reflection law leads to a concen- 
tration of the reflected energy density into a narrow ray 
tube upon reflection as displayed in Fig. G] Theoretical 
descriptions of this reflection process have been framed 
largely in terms of linear and stationary wave dynam- 
ics 0) H- However, when the slope angle, 7, is equal 
to the incident wave angle, (J, these restrictions lead to 
an unrealistic prediction: The reflected rays lie along the 
slope with an infinite amplitude and a vanishing group 
velocity. Theoretical results have recently healed this sin- 
gularity by taking into account the role of transience and 
nonlinearity 

Following preliminary oceanographic measurements |a, 
Q, Eriksen [3 has beautifully observed, near the bottom 
of a steep flank of a tall North Pacific Ocean seamount, an 
internal wave reflection process leading to a clear depar- 
ture from a Garett-Munk model 9] for wave frequencies 
at which ray and bottom slopes match. Several exper- 
imental facilities [H El 13 Hi 111 13 

have therefore 

been dedicated to the understanding of the internal wave 



FIG. 1: Schematic view of the reflection process when the 
incident wave nearly satisfies the critical condition 7 » /?. 
The group velocity of the reflected wave makes a very shallow 
angle with the slope. c 9 indicates the incident and reflected 
group velocities. The reflection law leads to a concentration 
of the energy density into a narrow ray tube. 
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reflection and associated instabilities. 



However, results of a controlled laboratory experiment 
close to the critical conditions are still lacking since pre- 
vious ones with a moving paddle GHHE! 

at one end 

of a very long tank, or by the vibration of the tank it- 
self 01 j does not generate a clear incident wave-beam as 
needed for a careful study. In addition, a direct compar- 
ison with the recent and complete nonlinear theory near 
the critical reflection would be possible. Finally, the goal 
is to improve the understanding of the possible mixing 
mechanism near the sloping topography of ocean as very 
recently initiated by MacPhee and Kunze 0] by exhibit- 
ing the instabilities mechanism leading to mixing. The 
paper is organized as follows. The experimental setup is 
described and carefully explained in Sec. [HI The main 
experimental results are presented in Sec. IIIII and com- 
parisons with the theory is also provided. Finally, Sec. II VI 
erieoiitiaMaf/ conclusions and perspectives. 
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II. EXPERIMENTAL SETUP 

The experimental setup consists of a 38 cm long Plex- 
iglas tank, 10 cm wide, filled up to 22 cm height with 
a linearly salt-stratified water obtained by the "double 
bucket" method The choice of salt, Sodium Nitrate 
(NaNOa) snow, is motivated due to its highly solubility 
in water leading to a salty water viscosity close to the 
fresh one. Moreover, this salt allows to reach a strong 
stratification: The fluid density ( 1 < p(z) < 1.2 g/cm 3 ) 
measured at different altitudes (22 > z > cm) by a 
conductimetric probe leads to a linear vertical density 
profile of slope dp/dz ~ —0.0104 g/cm 4 . A rectangular 
Plexiglas sheet, 3 mm thick and 9.6 cm wide (to allow 
exchange of water) is introduced at one end of the tank 
with an angle 7 = 35° with respect to the horizontal 
tank bottom (see Fig. [JJ to create the reflective sloping 
boundary. 

Internal waves are generated by a sinusoidal excita- 
tion provided by the vertical motion of a horizontal PVC 
plunging cylinder (3.1 cm in diameter and 9.4 cm long). 
The cylinder is located roughly midway between the base 
tank and the free surface. This wavemaker is driven by an 
electromagnetic vibration exciter powered by a low fre- 
quency power supply. Optical measurements confirm [l^ 
that the cylinder motion is sinusoidal without distortions 
for vibrational frequencies 0.2</<0.5Hz and maximal 
displacement amplitudes up to 8.5 mm (peak to peak). 

The well-known and non-intuitive dispersion relation 
of internal waves in an incompressible, inviscid and lin- 
early stratified fluid reads 

w = = ±A sin/3 , (1) 

I k I 

where k± is the component of k perpendicular to z- 
axis, N = ^/—(g/po)dp/dz is the constant buoyancy (or 
Brunt- Vaisala) frequency, p(z) the fluid density at alti- 
tude z, g = 981 cm/s 2 the acceleration of gravity, and 
po ~ 1 g/cm 3 a reference density. Thus, from Eq. QJ, 
the wave frequency, to = 2nf, determines the inclination 
angle (3 of the phase surfaces with the vertical, and not 
the magnitude of the wave vector k. From Eq. (JJJ , is 
also the angle between the group velocity c g = dcu/dk 
and the horizontal, since c g _L k. Thus, for a given strat- 
ification and frequency, internal waves of low (higher) 
frequency propagate at low (steeper) angle. 

The outward radiation of energy is thus along four 
beams oriented at an angle j3 with the horizontal, the fa- 
miliar St Andrews Cross structure [20j . The beam prop- 
agating directly toward the slope has been singled out by 
adding a grid on the surface of water. This grid strongly 
damps the three other wave beams that propagate to- 
wards the free surface, and thereby prevents reflection of 
such beams. A planar wave pattern consisting of parallel 
rays is thus generated from the wavemaker. Although it 
is well-known that the spatial spectrum of waves gener- 
ated by an oscillating cylinder is large [T5I l2lt |22T | , it has 



been experimentally checked that the dominant wave- 
length is approximately equal to the cylinder diameter 
(see below). Moreover, measurements confirm that the 
low vibrational amplitudes of cylinder do not influence 
strongly the wavelength generated. 



Different visualization methods are used to study the 
reflection mechanism of such internal gravity waves by 
a boundary layer. First, the usual shadowgraph tech- 
nique |23j allows visualization of the qualitative and 
global 2D evolution of the incident and reflected waves. It 
involves projecting a point source of light through strat- 
ified water onto a screen behind the tank. The optical 
refractive index variations induced by the fluid density 
variations, allows to observe isodensity lines (or isopyc- 
nals) on the screen, located perpendicularly to the light 
source and parallel to the longest wall tank. It is there- 
fore possible to measure the group velocity angle, /3, and 
the phase velocity v v of the incident wave. For various 
frequencies of excitation, 0.2 < f < 0.5 Hz, the angle f3 
of the St Andrews Cross is measured on a screen lead- 
ing to a linear relation between u> and sin/3 as predicted 
by Eq. JTJ) with a slope of iV = 3.1 ± 0.1 rad/s, for the 
stratified fluid prepared as above. This value is in good 
agreement with the above static one obtained from the 
density profile with a conductimetric probe. The cut- 
off frequency is then f c — N/(2%) ~ 0.5 Hz. By time 
of flight measurements between signals delivered by two 
photodiodes, 1 cm apart, each of 7 mm 2 area, located 
along the propagation direction, a value v v — 0.6 ± 0.4 
cm/s was obtained for an excitation frequency / = 0.25 
Hz. When these signals are cross-correlated by means of 
an spectrum analyzer, the averaged dephasing time leads 
to v v = 0.7 ± 0.3 cm/s, close to the previous value. The 
wavelength of the incident wave is thus A = v v / f ~ 3 
cm, corresponding as expected to the oscillating cylin- 
der's diameter. However, neither this shadowgraph tech- 
nique nor one p3 | using passive tracers (fluorescein dye) 
is sufficiently sensitive to observe quantitative and local 
internal wave properties closed to the reflective boundary 
layer. 

Accordingly, the classical Schlieren method [2^, [24| of 
visualization has been used. Let us just note that behind 
the tank, the light beam is refocused by a lens a small 
distance after a slit (instead of the usual knife blade to 
increase contrast) to filter the rays. The slit is oriented 
orthogonally to the slope generating straight horizontal 
fringe lines in the case of no excitation. The image of the 
observation field (strongly dependent of the 7.5 cm diam- 
eter lens) is focused on the screen by a last lens. The in- 
ternal wave, producing density disturbances, causes lines 
to distort, this distorting line pattern being recorded by 
a camera. Note that this experimental technique is sen- 
sitive to the index gradient, and therefore to the density 
gradient, orthogonal to the slit, i.e. parallel to the slope. 



3 



III. EXPERIMENTAL RESULTS AND 
COMPARISON WITH THEORY 

The Schlieren technique allows to carefully observe 
quantitative and local internal wave properties during the 
reflection process. Critical reflection arises when an inci- 
dent wave beam with an angle of propagation /3 reflects 
off the slope of angle 7 ~ (3, the reflected wave being then 
trapped along the plane slope. This corresponds to a crit- 
ical frequency f c = iVsin7 ~ 0.28 ± O.Of Hz, N being 
equal to 3.1 ±0.1 rad/s for all experiments. It is possible 
to observe that the isodensity lines (isopycnals), initially 
horizontal without excitation, are bent for an excitation 
near f c (0.78 < f/f c < 1.14), and fold over themselves 
along the length of the slope. Figures show a time 
sequence of constant density surfaces, depicting sequen- 
tial snapshots of the flow throughout one period of its 
development. These pictures show strikingly the distor- 
tion of the isopycnals in the slightly subcritical case with 
f / fc — 0.78. In panel (a), the density disturbance is very 
small and one can observe essentially the initial horizon- 
tal background stratification. This background stratifi- 
cation, usually obtained with dye fluorescein, should be 
invisible with this schlieren method. However, as pre- 
viously reported by MacPhee and Kunze [lfij. the com- 
parison between shadowgraph and Dye fluorescein visu- 
alizations allows to identify these lines with isopycnals. 
In panel (b), the disturbance generated by the incident 
wave breaking against the slope begins to 'fold-up' the 
isopycnals. As time progresses (see panel c), wave over- 
turning develops around a front: The buoyancy becomes 
statically unstable. This overturned region climbs along 
the slope as time continues, and the folded isopycnals col- 
lapse into turbulence that mixes the density field within 
the breaking region (see panel d). The maximum thick- 
ness of this reflected disturbance is of the order of 5 mm, 
and decreases with time as theoretically predicted |5| . Fi- 
nally, the flow begins to relaminarize (panel e). One can 
check that panels (a) and (f) are almost identical, show- 
ing that the flow is entirely restratified (panel f ) after one 
period of excitation. 

Figures|21have been analyzed with an image processing 
software (Scion-Image) to extract the isopycnals from the 
pictures. A typical experimental result for the distortion 
of isopycnals is reported in Fig. EJi, and is compared with 
a theoretical result in Fig. EJd. The analytic solution of 
the density field of the initial value problem in the critical 
case reads M 



P = Po 



N 2 I k| sin 2 (2/3 ) 
1 ZCOS7 — ipB\ 1 ■ — 

9 \ V 2w + 



(2) 



Ji being the Bessel function, lj + the positive solution 
of Eq. Q), ip the maximum amplitude of the stream- 
function and x the horizontal coordinate. Figure El show 
a good qualitative agreement between experimental and 
theoretical results, the value for the time t being arbitrar- 
ily chosen. Far from the slope, the density disturbance is 
very small and one sees essentially the initial background 
stratification. Closer to the slope, the disturbance folds 
up the isopycnals, and this produces a region of static 
instability. 

Recording several isopycnals and using image process- 
ing, it is also possible to follow the temporal evolution 
of a single isopycnal during its overturning. As this phe- 
nomenon is periodic with a period T = 1//, it is possible 
to reconstruct from this temporal evolution the density 
profile picture at a given time t. This allows to follow 
the position, and therefore the propagation velocity of 
the front along the slope at different times. The front is 
defined as the inflexion point (represented by the star in 
Fig. Ep) of the followed isopycnal. Figure^ shows, dur- 
ing two periods, the isopycnal front position along the 
slope as a function of time. The periodic evolution of 
this front position is clearly observed, and the local slope 
of curves in Fig. 0Ji allows to roughly measure the front 
velocity as a function of time, as reported in Fig. Ep. 

The front velocity from its creation to its collapsing in- 
creases from 0.5 cm/s up to 3 cm/s. The front has thus 
traveled roughly 4 cm in one period (~ 4.5 s). This leads 
to an averaged front velocity in agreement with the phase 
speed measurement obtained from the shadowgraph vi- 
sualizations. 

The wavemaker frequency is now increased up to / = 
0.32 Hz, to have an incident planewave tilted with an 
angle (3 steeper than the slope angle 7. In this slightly 
supercritical case (f / f c = 1.14), intrusions are still ob- 
served, but the density field doesn't fold up so abruptly 
and does not lead anymore to overturning instability (see 
Fig- E| . Except the value of the frequency /, all others 
parameter values have been kept identical to the ones 
in Fig. The instant of this snapshot has been chosen 
when the isopycnal distortion is the largest [25|. This 
distortion is clearly far from leading to an overturning in- 
stability as encountered in the subcritical case of Fig. [2;: 
The reflected wave is not trapped along the slope in the 
boundary layer, and consequently the isopycnals are not 
overturned. Moreover the density front velocity is mea- 
sured roughly constant, 0.5 cm/s, during two periods of 
vibration. Both differences confirm that the singularity 
appears only in the critical case 



IV. CONCLUSIONS AND PERSPECTIVES 



where 



B 



Ji sin (oj + t— I k I sin(/3 + 7)2;) , (3) 



J, = Ji (2 v/2w + cos 2 /3 k I tz 



(4) 



The Schlieren technique allows us to study the spa- 
tiotemporal evolution of the internal waves reflection 
close to the critical reflection, where nonlinear processes 
occur. The dynamics of isopycnals is then found in agree- 
ment with a recent nonlinear theory Moreover, this 
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(a) t = 



(b) t = T/5 




(e) t = 4T/5 



(f) t = T 




FIG. 2: Schlieren pictures showing the slightly subcritical reflection (f/fc = 0.78) of an internal wave on a slope, during one 
incident wave period T. The slope (thick black line) has an angle 7 = 35°. The incident wave plane comes in from the left 
(inclined black region near the top left corner between blue and yellow regions). The reflected wave plane is hardly noticeable. 
Wavemaker vibrational frequency and peak to peak amplitude are respectively / = 0.22 Hz and 6.7 mm. (Color pictures). 
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FIG. 3: (a) Experimental isopycnals extracted from a region 
of Fig.|2J:. (b) Theoretical isopycnals from Eq. The star 
indicates the position of the front. See text for parameters. 
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FIG. 4: Temporal evolution of the isopycnal front position 
(a) and velocity (b) along the slope, during two periods of 
vibration. T ~ 4.5 s, f/f c = 0.78 and A pp = 6.7 mm. 

experiment confirms the theoretically predicted scenario 
for the transition to boundary-layer turbulence responsi- 
ble for boundary mixing: The growth of a density pertur- 
bation produces a statically unstable density field which 
then overturns with small-scale fluctuations inside Il6l. 



Panel (d) of Fig.[2]is characteristic of the onset of turbu- 
lence triggered by overturning instability near the slope. 

This turbulent mechanism is likely responsible for the 
formation of highly "stepped" temperature profiles near 
steep slopes in lakes |26|. Moreover, the formation of sus- 
pended sediment layers, called nepheloid layers, at conti- 
nental slopes has been linked to critical angle reflection of 
internal waves |27J , and suggests that this reflection 
process plays an important role in the seawards trans- 
port of sediments in the fluid. A possible extension of 




FIG. 5: Isodensity lines during the slightly supercritical re- 
flection (f/fc = 1.14) of an internal wave on a slope. The 
slope (thick black line) has an angle 7 = 35°. The incident 
wave plane comes in from the left (inclined white region near 
the top left corner between blue and yellow regions). The 
reflected wave plane is hardly noticeable. Vibrational param- 
eters: / = 0.32 Hz and A pp — 6.7 mm. This picture should 
be compared with Fig. [3J:. (Color picture). 



this work, with a smaller slope angle to be closer to real 
oceanographic situations, would be to study and char- 
acterize the long time behavior and diffusion process of 
such particle layers toward the fluid interior. 
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